Logarithmic distribution of mean velocity and turbulent kinetic energy in a pipe flow 
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A Lie-group based similarity theory is developed for both momentum and energy distributions in a turbulent pipe 
flow, leading to asymptotic logarithmic profiles of mean velocity and turbulent kinetic energy. Both channel and pipe 
data over a wide range of Re yield 0.45 to be the universal Karman constant. A new spatial invariant characterizing 
outer dynamics is discovered and validated by reliable experimental data. The theory predicts the mean velocity profile 
(MVP) with 99% accuracy for high Re experimental data (up to 40 millions), and offers a quantitative explanation 
Q-t for recent observation of logarithmic kinetic energy distribution by Hullmak et al. (Phys. Rev. Lett. 108, 094501). 

Turbulent flows over objects form thin vorticity layers called boundary layers. As it is widely accepted that near-wall 
flow physics is autonomous and independent of the flow being external or internal, pipe flow forms an experimentally 
and numerically expedient canonical flow for the study of wall turbulence. Despite extensive efforts, the prediction of 
the mean velocity still relies on empirical functions 1] having limited accuracy and limited range of Reynolds numbers 
(Re). Hence, the problem continues to receive vivid attention with great experimental 0,01 and theoretical [1] efforts. 

From a statistical physics point of view, turbulent pipe flows are at a far-from-equilibrium state encompassing 
not only a cross-scale energy flux (cascade) but also momentum and energy fluxes in space. Understanding physical 
principles governing the non-homogeneous transport and non-uniform distribution of the mean momentum and kinetic- 
energy is a log-standing goal of the research. Nearly eighty years ago, Prandtl § and von Karman Q, independently 
proposed the concept of mixing length with a linear dependence on the distance from the wall, predicting a logarithmic 
MVP and hence friction coefficient. However, this empirical model has led to controversies: Barenblatt et al 0] have 
claimed that power-law is a better description; Goldenfeld _H has proposed a model for friction coefficient using 
•i-h ' a power-law description. A recent model of L'vov et al. [J] is particularly noteworthy, as its log-law description 
yields predictions of reasonable accuracy over a range of finite Re (see Fig. 2). Recently, a logarithmic scaling for the 
streamwise mean kinetic energy profile (MKP) is reported with no explanation. Clearly, a deductive theory for 
O r joint MVP and MKP is still missing. 

Here, we present a Lie-group based similarity theory for turbulent channel and pipe flows. The original idea was 
| presented in [l(| and formulated rigorously in llj. The goal of the theory is to find invariant solutions of the averaged 
flow equations based on a symmetry analysis of a set of new quantities, called order functions, which are introduced 
■^j- | in close analogy to order parameter [12j in the study of critical phenomena. Adding the order function to dependent 
. variables in the equations and then performing a dilation-group transformation yields a set of new, candidate invariant 
solutions, which defines a new method to measure the Karman constant. In this Letter, we compare the prediction 
to measured MVP in a turbulent pipe 13]. In addition, with a system similarity argument, we find a new spatial 
Q\ ■ invariant which predicts a logarithmic MKP at high Re. 

, Theory for mean velocity - In a pipe, the mean momentum equation (MME) is [i| : 

(N- 

S + + W + = — -(uv) + =t+ (1) 

•rH ■ 

where U + = (u) + is the streamwise mean velocity, S + = dU + /dy + is the viscous stress, W + = — (uv) + is the 
Reynolds stress, and = r is the total wall shear stress, r is the distance to the centreline, y + the distance to the 
wall, and u, v are streamwise and vertical fluctuating velocities, respectively, and superscript + denotes 'wall units' 
normalization with friction velocity and viscosity. We now search for group-invariant solutions of (jTJ) , by introducing 
the mixing length, t\ t = vW+/S + . A formal Lie-group analysis of JT]) adding £m (= £ M /Re T ) and its gradient, 
(Mm /dr, as new dependent variables [TTJ] shows three kinds of invariant solutions, of which the second kind corresponds 
to the situation of symmetry-breaking in l M but maintained symmetry in Mm I dr. This solution is expressed as: 

£ M (r)*-(l-r m ) (2) 
m 

where m is a scaling exponent characterizing the bulk flow, and k is the classical Karman constant, as can be seen 
when taking the limit to the wall, r — > 1, £m « ay, i.e. Karman's linearity assumption. 

Now, we establish a complete expression for the mixing length in the outer flow encompassing a bulk and a central 
core. Define another length involving S + , W + and dissipation e + , using the eddy viscosity v t = W + / S + [5J, which 
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FIG. 1: (color). Plot of measured U~[ versus theoretical f{r;r c ) (vertically shift by one) for five sets of data : DNS channel 
(black and red; r c = 0.27), DNS pipe (green; r c — 0.23), experiments pipe (solid diamonds, disks and squares; r c — 0.67). Note 
a good linear relation with slope 0.45 for all Res. The original plot is shown in the inset. 



gives £^ — (W + 1 S + ) 3 / 4 /e +1 / 4 = i^ 4 /e +1 / 4 . A comparison to the Kolmogorov dissipation length 77 = i/ 3 / 4 /e 1 / 4 
suggests that £ v signifies the energy input length scale, while 77 the output length scale of the energy cascade. Verify 
that 

£ M = H@ X '\ (3) 

where 9 = e + /(S + W + ) = e + /P + is the ratio between dissipation and production. is an important quantity, 
denoted to be an order function of the second kind [l(| • 

A two-layer model can be readily derived with the above definitions. In the bulk flow, the quasi-balance (lij 
corresponds to ~ 1. Near the center line as r — > 0, W + « r, S + « \fr/t M ~ r an< ^ e > 0, then Im ~ r -1 / 2 , 
but £ v — >• C = 0(1) and ~ r~ 2 . This transition is physically due to the switch of the generating mechanism of 
fluctuations - from mean shear production to turbulent transport (vanishing mean shear at the center). A simple 
ansatz taking into account of the transition of from r~ 2 to 1 is = 1 — c + cr~ 2 with a constant c, or 

= 1±M^ (4) 
1 + r 2 ' {) 



where the parameter r c — yc/(T— c) represents a critical radius characterizing the core region. In this core region 
(r c > r —> 0), l v ~ £q and w r _2 r 2 /(l + r 2 ) together determine the central behavior of Im from ©. 

In the bulk flow (r > r c ), ~ 1, £^ ~ A phenomenology yields an estimation of m in @ as following. 

Consider the scaling property of l v . Similar to Im, only dl v /dr has a scale-invariance property: d£ u /dr oc r" l_1 , 
which corresponds to the Lie-group similarity of the second kind [ll|. The exponent m can be derived by assuming 
that: d£ v j dr oc J P + rdr, i.e. the volume integral of the turbulence production, which equals the amount of kinetic 
energy converted from the mean flow. The validity of this assumption relies on an intriguing connection between the 
mean flow and fluctuation energy yet to be uncovered, which we defer to future study. Since P + = S + W + oc r 2 and 
using the wall condition (l v = as r — > 1), we obtain l v = £o(l — r 5 ). Furthermore, taking the limit r — > 1 and using 
Im ~ ny and ~ 1, we obtain an important relation: k — 5£o- Hence, £ v — k(1 — r 5 )/5. This derivation immediately 
predicts that, for channel flow, m = 4, since dl v jdr oc J P + dr with a flat plate. Substitute the expression of £ v and 
(@| into (|3]), we thus obtain: 

£M(r) = ^-{l-r m )(l + (r c /rf) 1/i (5) 

where Z c = (1 + r 2 ) 1 / 4 , and m = 5 (4) for pipe (channel) flows. 

Measurement of k and prediction of the MVP - In the outer flow, W + w £^S +2 w r, which yields S + w \/t/£ m . 
Integrating it using ([5]) yields an expression for the mean velocity defect, i.e. — — U + (r) = Re T S + dr = 
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(l/«)/(r; r c ), where is the mean velocity at the centreline, and 



r < 



f(r; r c ) = mZ c ^ (i - - = ^ . (6) 

The linear relation between f(r;r c 

) and Uj can be subjected to experimental tests, with a fixed r c . In FiglTJ 
theoretical f(r;r c ) versus empirical (measured) Uf is plotted for a wide range of data, from direct numerical sim- 
ulations (DNS) of channel and pipe flow at moderate Re fl5l - fl7j and experiments of pipe flows [l3| at very high 
Re. The linearity is remarkably observed with a slope of 0.45 ± 0.01, consistent with our earlier study of chan- 
nel flows jl8j • Note that r c can be determined in a rational way, by an evaluation of the minimum relative error, 

°v = jf J2 (l — K Ud EXP ( r i)/ f( r i'i r c)) , where k is the value from the least-square fitting above. In practice, the 
relative error is evaluated in a domain defined as 200 < y + < 0.9Rc r . This procedure yields a r c pa 2/3 for Princeton 
data at high Re T > 5000, while at the moderate to low Re, the DNS data [H-[l3] show a smaller r c pa 0.23 ~ 0.27 
for Re T < 1000 (both channel and pipe). With r c so determined, k (i.e. the slope) is measured with high confidence: 
all data show that K pa 0.45. 

This value of k is 10% higher than generally accepted value (0.41) [l[. Initially, it was a bit surprising, but a careful 
scrutiny confirms the internal consistency of the procedure. Note that in our definition, n is a coefficient defining the 
outer flow, and its measurement involves little ambiguity. The fact that the measured value is accurately constant 
for both channel and pipe flows and for a wide range of Re's, and that it is exactly the historic Karman constant in 
the overlap region, suggests that (J5J) is a better definition for n in channel and pipe flows, and taking into account 
ex plic itly the form of the outer flow makes the measurement of K more robust, compared to previous measurements 
(3. Il9|. Another interesting prediction is an asymptotic centreline dissipation at high Re: 

4 P iMCH) = lim{s+w+e) = „ ±|M (7) 

r^o y ' K Z c 3 Re T Re r v ' 

The constant 3.7 for channel flow becomes 1.18 at Re T pa 940 (with r c ~ 0.27), which agrees with measured value 
e+ DNS pa 1.05 /Re T from DNS data of [H]. 

In order to predict the MVP, one needs an additional constant. Using ([6]), we can express the MVP as 

U + (r) = U+-f(r;r c )/K (8) 

Our analysis of the Princeton pipe data at high Re show that rj lni?e T /0.45 + B c , with B c pa 8.3. With three 
parameters: n ss 0.45, r c pa 0.67 and B c pa 8.3, the high-i?e outer flow MVP is completely specified by (|6|) and ([8]). In 
addition, an asymptotic calculation of ([S]) gives lim /(r; r c ) = lim pdJ~f (r) pa — In (1 — r) + kB^ with B& ?a 1.7, which 

i — 7-1 r— ^1 

only depends on r c . Then, ((BJ) yields an approximate log-law in the overlap region: 

U + (r) sa lny+/0.45 + 6.6 (9) 

where the additive constant is found from B c — B^ pa 8.3 — 1.7. In Fig(21 the theoretical MVP © is shown to agree 
with the Princeton pipe data for the entire profile with 99% accuracy for Re up to 4 x 10 7 (better than [J]). This 
unprecedented accuracy settles the debate between the logarithmic law and power law with a rational description of 
the bulk flow; it further shows that turbulence in pipe flows indeed admits an analytic solution. 

Predictions for fluctuations - The mean kinetic-energy equation (MKE) [14| can be rearranged in a form similar to 
©as: 

- ^pr + < um; > + = t k (io) 

where K + = (uu) + is the streamwise kinetic energy, and = J y (S + W + — e^)dy + — (pu) + involves the integration 
of turbulent production, SW, dissipation, e#, and a term due to pressure transport, (pu) + . Note that one often 
assumes a quasi-balance, then, would be small. Detailed examination of empirical data shows that the integrated 
deviation from the quasi-balance, although small, is fully responsible for the non-uniform distribution of the kinetic 
energy. 

The similarity between © and (|10[) is a major focus of this paper. Denote = —dK + /dy + and = (uuv} + , 
thus, the MKE has a similar form as 



(11) 
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FIG. 2: (color). Theoretical (solid lines) and measured MVPs, which are staggered vertically by five units for clarity. Inset 
shows the relative errors, (U EXP /U THE — 1) x 100 % - our predictions (red solid symbols) are uniformly within 1%. Also 
included (blue open symbols) is the model of L'vov et al. 0], which reveals systematic deviations at high Re due to inappropriate 
boundary constraints [20l ]. 



The two terms on the l.h.s. represent the viscous diffusion and turbulent transport of K . Analysis of DNS data of 
channel flows [l6| reveals indeed two similarities (results not shown): first, all terms go to zero near the center line, 
and W + ~ W£ and S + ~ Sp, second, W+ > S%, but W + /S + ~ W£/S^. These similarities suggest that © 
and pop may differ by a constant factor a. Assuming this is true, multiplying ([1]) by a and subtracting (|10|) yields 
d{aU+ + K+)/dy+ « 0, leading to 

aU+ (r) + K+ (r) w aU+ + K+ C (12) 
where K J is central kinetic energy, with 

C = (a/«) lni?e r + aB c + K+ (13) 

being a constant in the outer flow. With experimentally measured U + (r) and K + (r) , the validity of (fT2"j) is successfully 
tested (Fig[3]) with a a « 0.56 (Fig[5|) by a linear fitting at small r (where the linearity is accurate because both terms 
have a quadratic dependence on r) . Figf5] shows a clear evidence of a spatial invariant over an increasing radial 
domain with increasing Re: the extent reaches almost the entire radius at high Re. In addition, the constant C can 
be measured directly from data at each Re, shown in the inset: C ~ 1.251ni?e T + 5.4. This measurement yields 
K+ w 0.7. 

A specific prediction of (fl"2"]) is that at high Re, K must have a logarithmic profile (with a negative sign), since U 
has a logarithmic profile. In the overlap region, K + w — (a/n) my + cuBd + Kj w —1.25 In y + 1.65, which reproduces 
well the empirical observation of Hulkmark et al. [9]: K + sw —1.25 In y + 1.61. Finally, we predict the outer profile of 
the MKP as 

K+ & a f(r;r c )/K + K+. (14) 

As shown in Fig(H it agrees well with empirical data, especially at high Re, better than that proposed by Alfredsson 
etal[2l|. 

In summary, we have achieved a simultaneous description of the outer MVP and MKP by exploring the similarity in 
the MME and MKE, which discovers a new spatial invariant in the radius direction of a turbulent pipe. The Lie-group 
based theory yields a new procedure for measuring k and a universal value of 0.45 for both channel and pipe flows. 
The predicted MVPs achieve a 99% accuracy compared to Princeton pipe data for a wide range of Re's. Note that 
the analysis has been successfully extended to incompressible, compressible and rough-wall turbulent boundary layers 
[22l [23l ] , and turbulent Rayleigh-Benard convection (temperature) . Those results will be communicated soon. 

We thank Y. Wu for helpful discussions. This work is supported by National Nature Science Fund 90716008 and 
by MOST 973 project 2009CB724100. 
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FIG. 3: (color). Evidence of linear dependence between U~f and K + . The coefficient a ~ 0.56 is obtained from fitting the 
largest Re data. The inset validates the invariance law (If 21) using experimental K [9J and theoretical U at the same Re with a. 
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FIG. 4: (color). Comparison between the theory (|f 4[l (solid lines) and experimental data (symbols) [|| of the MKP. 
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